Localization Transition of the 3D Lorentz Model and Continuum Percolation 
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The localization transition and the critical properties of the Lorentz model in three dimensions are investi- 
gated by computer simulations. We give a coherent and quantitative explanation of the dynamics in terms of 
continuum percolation theory and obtain an excellent matching of the critical density and exponents. Within a 
dynamic scaling Ansatz incorporating two divergent length scales we achieve data collapse for the mean-square 
displacements and identify the leading corrections to scaling. We provide evidence for a divergent non-Gaussian 
parameter close to the transition. 
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Transport in heterogeneous and disordered media has im- 
portant applications in many fields of science including com- 
posite materials, rheology, polymer and colloidal science, and 
biophysics. Recently, dynamic heterogeneities and grow- 
ing cooperative length scales in structural glasses have at- 
tracted considerable interest |lL|2|]. The physics of gelation, 
in particular of colloidal particles with short range attraction 
is often accompanied by the presence of a fractal 
cluster generating sub-diffusive dynamics. It is of fundamen- 
tal interest to demonstrate the relevance of such heterogeneous 
environments on slow anomalous transport. 

The minimal model for transport of particles through a ran- 
dom medium of fixed obstacles, is known as Lorentz model, 
and already incorporates the generic ingredients for slow 
anomalous transport. Earlier, the Lorentz model has played 
a significant role as a testing ground for elaborate kinetic the- 
ories, shortly after the discovery of long-time tails in auto- 
correlation functions for simple liquids in the late 1960s |7j, 
since the non-analytic dependence of transport coefficients on 
frequency, wavenumber, and density predicted for simple liq- 
uids l8|l9l flolfTllfl2il has a close analog in the Lorentz model 

Bm. 

The simplest variant of the Lorentz model consists of a 
structureless test particle moving according to Newton's laws 
in a (i-dimensional array of identical obstacles. The latter are 
distributed randomly and independently in space and inter- 
act with the test particle via a hard-sphere repulsion. Conse- 
quently, the test particle explores a disordered environment of 
possibly overlapping regions of excluded volume; see Fig.^ 
Due to the hard-core repulsion, the magnitude of the particle 
velocity, v = |v|, is conserved. Then, the only control param- 
eter is the dimensionless obstacle density, n* := na d , where 
a denotes the radius of the hard-core potential. At high densi- 
ties, the model exhibits a localization transition, i. e., above a 
critical density, the particle is always trapped by the obstacles. 

Significant insight into the dynamic properties of the 
Lorentz model has been achieved by a low-density expansion 
for the diffusion coefficient by Weijland and van Leeuwen 
fklll rigorously demonstrating the non-analytic dependence on 
n*. As expected, for low densities the theoretical results com- 
pare well with Molecular Dynamics (MD) simulations fT3l . 



Elaborate self-consistent kinetic theories flfl \vfa have al- 
lowed going much beyond such perturbative approaches. 
They give a mathematically consistent description of the lo- 
calization transition, which allows to calculate the critical den- 
sity within a 20% accuracy and extend the regime of quanti- 
tative agreement to intermediate densities. In addition, they 
have provided a microscopic approach towards anomalous 
transport and mean-field-like scaling behavior flr^l . 

A different line of approach focusing on the localization 
transition starts from the fractal nature of the void space be- 
tween the overlapping spheres in the Lorentz model and con- 
siders it as a continuum percolation problem fisl fl^. Eol EH 
12211 . which in this context has also been termed "Swiss cheese" 
model 1 20] . These authors conjectured that the transport prop- 
erties close to the percolation threshold can be obtained by 
analyzing an equivalent random resistor network. The equiv- 
alence, however, has been shown only for geometric proper- 
ties close to the percolation point |23]. As a peculiarity of 
continuum percolation, differences to lattice percolation may 
arise due to power law tails in the probability distribution of 




FIG. 1: (color) Typical particle trajectories in a 2D Lorentz model 
slightly below n* over a few thousand collisions each. Colors en- 
code different initial conditions; obstacles have been omitted for clar- 
ity. Most trajectories being in the percolating void space have some 
overlap; a few trajectories are confined to finite clusters. Blow-up: a 
particle squeezes through narrow gaps formed by the obstacles. 
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FIG. 2a: (color) MSD 5r 2 {t) for various obstacle densities n* vary- 
ing from 0.30 (top) to 1.10 (bottom). The thick black line represents 
a power law, 8r' 2 (t) ~ t 2 ^ z with z = 6.25. 
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FIG. 2b: (color) Scaling functions 5f±(l) for the MSD. Right panels 
include corrections to scaling at leading order. Units are chosen such 
that v = a = 1; color key as in Fig.l2al 



the conductances ("narrow gaps"). Such random resistor net- 
works have been investigated extensively by means of Monte- 
Carlo simulations 1241 12511 and renormalization group tech- 
niques |26, 27], providing reliable numeric and analytic re- 
sults for the critical behavior [28]. 

In this Letter, we present a direct numerical analysis of 
the dynamic properties of the Lorentz model without resort- 
ing to random resistor networks. By means of extensive MD 
simulations, we obtain a quantitative description of the dy- 
namic properties over the full density range, in particular, 
focusing on both sides of the critical region. This allows 
for a quantitative test of the conjectured mappings to contin- 
uum percolation theory. Furthermore, we explore the range 
of validity of the dynamic scaling hypothesis for the Lorentz 
model [29]. The probability distribution of particle displace- 
ments, i. e., the van Hove self-correlation function, G(r, t) := 
\S(r— AR(t))), and its second moment, the mean-square dis- 
placement (MSD), Sr 2 (t) := (|AR(t)| 2 ), are the appropriate 
quantities for this purpose; AR(i) = R(i) — R(0) denotes 
the displacement of the test particle at time t. 

Over a wide range of obstacle densities, we have simulated 
several hundred trajectories in three dimensions, employing 
an event-oriented MD algorithm. For each of N r different 
realizations of the obstacle disorder, a set of N t trajectories 
with different initial conditions is simulated. Below the crit- 
ical density, we have chosen N r > 25 and Nt > 4. At very 
high densities, where the phase space is highly decomposed, 
these values have been increased up to N r x N t — 600. In 
order to minimize finite-size effects, the size of the simula- 
tion box, Lbox, has been chosen significantly larger than the 
correlation length £, L box = 200a > £ Q. 

The results for the MSD cover a non-trivial time window 
of more than seven decades for densities close to the transi- 
tion, see Fig.|2a| At low densities, one observes only a trivial 
cross-over from ballistic motion, Sr 2 (t) = v 2 t 2 , to diffusion, 
5r 2 (t) ~ t, near the mean collision time r = l/nnva 2 as 



expected from Boltzmann theory. With increasing density, 
an intermediate time window opens where motion becomes 
sub-diffusive, 5r 2 (t) ~ t 2 / z with z > 2. This time win- 
dow extends to larger and larger times upon approaching a 
certain critical density n*. For the density n* — 0.84, the 
sub-diffusive behavior is obeyed over more than five decades 
and is compatible with a value of z w 6.25. The power law, 
5r 2 {t) ~ t 2 / z '", indicated in Fig. [2a] discriminates nicely tra- 
jectories above and below n*. One also observes a density- 
dependent length scale £ characterizing the end of the sub- 
diffusive regime by Sr 2 (t) ~ £ 2 ; upon approaching n* this 
cross-over length scale £ is found to diverge. For long times, 
the dynamics eventually becomes either diffusive or localized 
for densities below or above n*, respectively. 

The diffusion coefficient D has been extracted from the 
long-time limit of 5r 2 (t)/6t; in Fig. [3] D is shown in units 
of the Boltzmann result, Do = tv 2 /3. With increasing den- 
sity, D is more and more suppressed until it vanishes at n* c as 
a power law, D ~ \e\^, where e := (n* — n* c )/n* c defines 
the separation parameter. Anticipating the exponent /i from 
percolation theory, a fit to our data yields the critical density, 
n* = 0.839(4) 113 511 . and the power law behavior is confirmed 
over five decades in D. Above the critical density, the long- 
time limit of the MSD is compatible with a power law over 
more than one decade, £ ~ e~ u+l3 / 2 , where v — (3/2 rj 0.68 
(bottom inset in Fig. [3}- Our finding of n* coincides with the 
percolation point of the void space 1 1 Si Il9l 13011 . This pro- 
vides clear evidence for the intimate connection between con- 
tinuum percolation and the Lorentz model, i. e., diffusion is 
not blocked as long as there is an infinite path through the 
medium — a purely geometric reason. 

Considering the underlying continuum percolation prob- 
lem, a geometric transition occurs at n*, above which the void 
space falls completely apart into finite clusters. Just below this 
density, the volume fraction P of the percolating void space 
(infinite cluster) vanishes as a power law, P ~ |e|' 3 . There are 
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FIG. 3: Suppression of the diffusion coefficient D/Do with increas- 
ing density n* . Top inset: Power-law behavior of D close to n* . Bot- 
tom inset: Divergence of the localization length t upon approaching 
n* with exponent v — (3/2 = 0.68. 



two divergent length scales characterizing the structure of the 
percolation network: the linear dimension of the largest finite 
clusters, £ ~ | sr | *^ , and the mean cluster radius (radius of gy- 
ration), The geometric exponents (3 and 
v are believed to be the same for lattice and continuum perco- 
lation 1 19]. Our results in Fig.[3]clearly identify the geometric 
mean cluster radius £ with the localization length of the MSD 
as anticipated by our choice of notation. 

In continuum percolation, transport of a particle is limited 
by narrow gaps in the void space. It was argued that this 
feature of the dynamics is captured by an associated random 
resistor network with a distribution p(W) of weak conduc- 
tances W exhibiting a power-law tail, p(W) ~ W~ a , a < 1 
for small W 12(1 12111 . Depending on the value of a, the 
suppression of diffusion, D ~ |e| M , may be dominated by 
this tail, and dynamic exponents become different from lat- 
tice percolation, fi > /i lat . In this case, the hyperscaling re- 
lation, fx = (d - 2)v + 1/(1 - a), holds ^ Q. There 
is a discrepancy in the literature about the value of a in the 
Lorentz model 12(1 l2ll 12811 . Only the result of Machta and 
Moore liUl . a = (d — 2)/(d — 1), is consistent with our 
data. In d — 3, it implies p, — v + 2 « 2.88, and there- 
fore, /i > /.t lat w 2.0 j3al . By means of a scaling relation l3lll . 
z = (2i/ — (3 +fj,)/(u — /3/2), one finds the dynamic exponent, 
z « 6.25, describing anomalous transport at critical density, 
8r 2 (t) ~ t 2 / z . Note that this dynamic exponent is not inde- 
pendent but entirely determined by the geometric properties 
of the random environment. 

In conclusion, the values obtained from the simulated MSD 
for the critical density n* , the dynamic exponent z as well as 
the exponents for the diffusion coefficient p and the localiza- 
tion length v — (3/2 agree with the predicted values for contin- 
uum percolation. Within the statistical accuracy, no deviations 
can be inferred. 

The quality of our data allows to go beyond determining 



critical exponents and to give a full analysis of the dynamic 
scaling properties. It has been argued by Kertesz and Met- 
zger l29ll that the van Hove correlation function obeys scaling. 
Rewriting their Ansatz in a more transparent way yields, 

G(r,t;s)=C 0/ "- d g±{r/Z,te-% (1) 

where Q± are master functions above (+) and below (— ) the 
critical density. This Ansatz clearly reflects the role of the 
two length scales: the correlation length £ rescales geome- 
try whereas the cross-over length scale £ rescales time. The 
scaling form of the MSD is easily inferred from Sr 2 (t; e) = 
Jd d rr 2 G(r,t;e) as, 5r 2 (t;e) = t 2 / z Sf 2 ± (t), where i ~ 
U~ z . Plotting t- 2 / z 8r 2 {t-e) versus t for various densities 
(left panels of Fig. I2b> . the data collapse nicely in the dif- 
fusive and localized regimes (t ~S> 1) and converge rapidly 
to the corresponding large-t asymptotes, (5f 2 _(t) ~ £ 1_2 / z 
and Sr+(t) ~ t~ 2 l z . Convergence to the critical asymptote, 
Sfj_ (t) ~ const, for t«l becomes increasingly better as the 
critical point is approached. 

Deviations from scaling can be rationalized by consider- 
ing the again universal corrections to scaling. Extending 
the Ansatz, Eq. Q, by an irrelevant parameter u leads to 
Sr 2 (t;e,u) = t 2 / z <&±{t£~ z , ut~ y ), where y is a universal 
exponent. Since $± is assumed to be analytic for small ar- 
guments, one obtains the leading order correction upon ex- 
panding <1>± to first order in u, 

5r 2 (t;e) = t 2 ^ z 6f 2 ± (t)(l + A ± (t)) , (2) 

introducing some analytic functions A±(t). Specializing 
Eq. (I3 to the critical density, i.e., t = 0, yields Sr 2 (t;e — 
0) oc t 2 / z {\ + Ct-y), with a single amplitude C = A±(t = 
0); it also identifies y as the leading non-analytic correction 
exponent at criticality. Our data for ri* = 0.84 are compatible 
with values for y between 0.15 and 0.4. For the following, 
we found the choice y = 0.34 and C = —0.8 reasonable, the 
value for y is supported by theoretical arguments to be pre- 
sented elsewhere. 

Inspection of Fig. l2bl re veals that corrections to scaling are 
less relevant for long times, t > 1, whereas significant devi- 
ations are visible in the critical regime, t 1. This obser- 
vation is consistent with the scaling behavior of the diffusion 
coefficient and the localization length, see Fig. [3] These find- 
ings also suggest approximating the corrections by its value 
at t = 0, i.e., substituting A±(t) = C in Eq. (0 for all 
times, 5r 2 (t;e) = t 2 / z 8f 2 ± {t) (l + Ct-y). With y and C 
already inferred from the data close to criticality, the correc- 
tion terms should apply for all densities. Indeed, including 
this leading-order correction improves the data collapse sub- 
stantially (Fig.|2bJ right panels). 

The presence of two different length scales, £ and £, in 
the scaling hypothesis, Eq. Q, is not manifested in the 
MSD; it will, however, affect the higher moments of the 
probability distribution, e. g., the mean-quartic displacement 
(MQD), Sr 4 (t; e) = Jd d r r 4 G(r, t; e). Above n* c , it is eas- 
ily inferred that the long-time limit of the MQD scales as, 
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FIG. 4: (color online) MQD 5r 4 (t) for densities above and below 
n*. The dashed and dotted lines compare the two exponents z and z. 
Inset: NGP a?2 it) for the same densities, its long-time limit increases 
as n* is approached. 

8r A (t — > oo) ~ £ 2 £ 2 . At the critical density, we obtain 
the long-time asymptote, 5r 4 (t) ~ t 4 / z , with the exponent 
~z := (2v - (3 + n)/{v - 13/ A) w 5.45 different from z. We 
have evaluated the MQD by our simulation and find agree- 
ment with the prediction of continuum percolation at a similar 
level of significance as for the MSD, see Fig.|4] In particular, 
for the density n* — 0.84 the MQD follows a power-law with 
the predicted exponent z for a time window of more than four 
decades. 

A more sensitive quantity is the (first) non-Gaussian pa- 
rameter (NGP), a 2 (i) := |<5r 4 (t)((5r 2 (i))" 2 - 1, quantify- 
ing deviations from a Gaussian distribution 13 311 . At critical- 
ity, it diverges as 0:2 (£) ~ i 4 / z ~ 4 / 2 ps ^0.097. direct obser- 
vation of this very small exponent is expected to be a con- 
siderably difficult task. The long-time limits of ct2(t) di- 
verge upon approaching n* from either above or below as 
otzit — » 00) ~ |e| - ^. In particular, the NGP does not van- 
ish in the diffusive regime close to the transition due to the 
presence of localized particles even below n'* . Although there 
are significant statistical errors in the data for the NGP, Fig.|4] 
(inset) provides evidence for a significant increase of ct2(t) as 
density approaches n* from either side. The properties of the 
NGP demonstrate that the presence of two divergent length 
scales is crucial for the understanding of the dynamics close 
to the localization transition. 

We are indebted to W. Gotze for valuable discussions and 
critical comments as well as W. Kob and A. Moreno for help- 
ful correspondence. 
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